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ABSTRACT 

The single-particle densitity of states (DOS) for the Pauli and the Schrodinger Hamil- 
tonians in the presence of an Aharonov-Bohm potential is calculated for different values 
of the particle magnetic moment. The DOS is a symmetric and periodic function of the 
flux. The Krein-Friedel formula can be applied to this long-ranged potential when reg- 
ularized with the zeta function. We have found that whenever a bound state is present 
in the spectrum it is always accompanied by a resonance. The shape of the resonance is 
not of the Breit-Wigner type. The differential scattering cross section is asymmetric if a 
bound state is present and gives rise to the Hall effect. As an application, propagation of 
electrons in a dilute vortex limit is considered and the Hall resistivity is calculated. 
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1 Introduction 



In this letter, nonrelativistic physics described by the Schrodinger and the Pauli equations 
is considered in the presence of an Aharonov-Bohm (AB) potential A(r) [|TJ. We shall use 
the regular radial gauge, in which 

$ a 

A r = 0, Ay — - — = - — $ - (1) 
Z7cr 2irr 

Usually, $ = a $o is the total flux through the flux tube and a > is the total flux $ in the 
units of the flux quantum $ , $o = hc/\e\. However, the AB potential can be considered 
in a more general sense, since, formally, the same potential (of nonmagnetic origin) is 
generated around a cosmic string. The parameter $ is then 1/Qniggs, ot = e/Qniggs, and 
$ = 2n/e (in the units h = c = 1) with e and Quiggs being respectively the charge of a 
test particle and the charge of the Higgs particle 0. In what follows a will be written as 
a = n + rj, where n = [a] is the nearest integer smaller than or equal to a and rj being the 
fractional part. The case of a nonsingular flux tube of finite radius R will be discussed, 
too, as it is important from the experimental point of view. Indeed, flux tubes realized in 
experiments such as vortices in superconductor of type II are never of zero radius. Our 
main results are: 

• the validity of the Krein-Friedel formula || £§] for the density of states (DOS) is for 
the first time established for a singular potential and the change Ap a (E) over all 
space of the DOS induced by the AB potential is calculated; 

• a resonance is predicted to occur whenever a bound state is present in the spectrum; 

• in contrast to zero modes ||, the number of bound states does depend on the 
regularization of the interior of a flux tube; 

• in the presence of a bound state, the differential scattering cross section is asym- 
metric and the Hall effect occurs. 

Details of our calculations and complete proofs are given elsewhere || |7], [UJ . Here, main 
ideas are presented and some of the proofs are outlined. 
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The Krein-Friedel formula (|T6|) gives the DOS as the sum over phase shifts and thereby 
relates the DOS directly to the scattering properties. It is therefore very useful to have 
its extension in the direction of singular (especially Coulomb) potentials. Here the Krein- 
Friedel formula is applied directly and consistency with previous results is shown. For ex- 
ample, in the particular case without bound states we confirm the anticipation of Comtet, 
Georgelin, and Ouvry || that the change of the DOS is concentrated at zero energy. 
The DOS provides an important link between different physical quantities: the partition 
function, virial coefficients, effective action, and in the relativistic case between induced 
fermion number and anomaly || lQfl . We have already used its knowledge to calculate 
the persistent current of free electrons induced in the plane by the AB potential (Tl| . 



The discovery of a resonance was quite unexpected. Rather suprisingly, the shape 
of the resonance [can be read from Eq. fl2"2])] is not of the Breit-Wigner form. Since 
the latter is a direct consequence of analyticity, it poses an interesting question on the 
analytic structure of scattering amplitudes for singular potentials. The resonance can 
have a profound influence on the transport properties of electrons in an experimental set- 
up where electrons can penetrate the interior of the flux tube provided that the latter is 
prepared is such a way that its interior is not isolated from the system under consideration. 
A realistic physical realization of the penetrable flux tube is that suggested originally by 



Rammer and Shelankov [12] and later realized experimentally by Bending, Klitzing, and 



Ploog jnj, i. e., to put a type II superconducting gate on top of the heterostructure 
containing the two-dimensional electron gas (2DEG) (see Fig. [I]). When a magnetic field 
is switched on the conventional superconductor is penetrated by vortices of flux with 
a = 1/2. Therefore, electrons from the heterostructure do not move in the homogeneous 
magnetic field but in the field of a penetrable flux tube. 

By considering the case of a regular flux tube with a finite radius R we shall show that 
a bound state can only occur if the gyromagnetic ratio g m is anomalous and greater than 
two. If g m equals exactly to two, then according to Aharonov-Casher theorem || zero 
modes occur. If g m < 2, the coupling with magnetic field is not sufficiently strong enough 
to form neither zero modes nor bound states. In the region g m > 2, i.e., exactly where the 
gyromagnetic ratio of electron (g m = 2.00232) lies, the coupling with the magnetic field 
is enhanced and zero modes turn out to be bound states. However, the number of bound 
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Figure 1: Black layer is a superconductor of type II put on top of the heterostructure 
containing the two-dimensional electron gas (dotted region). When this sample is put in 
a homogeneous magnetic field, magnetic field penetrates the superconductor in Abrikosov 
vortices. Therefore, electrons from the heterostructure do not move in the homogeneous 
magnetic field but in the field of a (penetrable) flux tube. 

state is generally higher than the number of zero modes and the number of bound states 
does not only depend on the total flux but also on the energy of the magnetic field. 

The differential scattering cross section is a periodic function of the flux a and asym- 
metric with respect to ip — > —ip, where p is the scattering angle J7|. The asymmetry of the 
differential scattering cross section is easy to understand as for a > bound states occur 
only for I < 0, I being the orbital angular momentum. The asymmetry of the differential 
scattering cross section has direct experimental consequences since it leads to the Hall 
effect. The Hall resistivity is calculated in the dilute vortex limit. 

One has the unitary equivalence between a spin 1/2 charged particle in a 2D magnetic 
field and a spin 1/2 neutral particle with an anomalous magnetic moment in a 2D electric 



field fbl and our results apply to this case as well. 



2 The Pauli and the Schrodinger Hamiltonian in the 
Aharonov-Bohm potential and scattering phase shifts 

Let us consider the Pauli Hamiltonian, 

(p _ jia) 2 
2m 

where fi = fis/s is the magnetic moment operator, s is the spin operator, and s is the 
magnitude of the particle spin. For electron \i e = —g m \e\%/4mc = —fiB9m/2, being 
the Bohr magneton, and g m is the gyromagnetic ratio that characterizes the strength 
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of the magnetic moment [15]. By separating the variables, assuming e 



H is 



written direct sum, H 
space L 2 [(0, oo), rdr] ]1|, [16 



1H1, of channel radial Hamiltonians Hi in the Hilbert 



d 2 



dr 2 



Id v a . . 

-~T + -o + 9m-s z 0[r). 
r dr r z r 



(3) 



Here v — |Z + a| and s z = ±1 is the projection of the spin on the direction of the flux 
tube H], jlBf . The Schrodinger equation is recovered upon setting s z = 0. For positive 
(negative) energies the eigenvalue equation in the /-th channel reduces to the (modified) 
Bessel equation of the order v = \l + a\, 



Hi ipi = k 2 ipi 



(4) 



with k = ^/2mE/h. In case of the impenetrable flux tube the spectra of both the Pauli 
and the Schrodinger equations are identical. There are neither zero modes nor bound 
states in this case [[jj . The boundary condition selects only regular solutions at the origin 
and the "spectrum" is given by 



ipi(r,<p) = J\i +a \{kr)e 



Hip 



Phase shifts 0, 



5, 



-7T 



\l + a\), 



(5) 



(6) 



are in general singular: they do not decay to zero in the limit E — > oo. 

In general case, Hamiltonians Hi for which \l + a\ < 1 admit a one-parametric family 
of self-adjoint extensions |H| [17 ]. They correspond to different physics inside the flux 
tube. The situation will be considered when bound states 



B t (r,(p) = K\ l+al (Kir)e 



1 1 ^ 



(7) 



of energy E\ = —{h 2 /2m)K 2 are present in the / = — n, —n — 1 channels, with n = [a]. 
In the presence of the bound states the scattering states (|D in these channels have to be 
modified. They become 



?pi(r,Lp) = J\i +a \{kr 



AiJ_\i +a \{kr) 



(8) 
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This is because Hi has necessarily to be a symmetric operator what already determines 
Ai to be 

M = (k/ Kl ) 2 », (9) 

i.e., energy dependent 0. The radial part of the general solution (^) behaves for (r — > oo) 
as 

(1 _ A, P in \ l + a \ \ 
l-A ie -'^ l+a \ J v ; 

Therefore, 

1 /i7i ,, i \ / sin(|i + ak) \ 

* = o< l ~ I + « + arctan " 1 — I , (11) 

2 \cos(|/ + a\7T) — A l J 

which determines S matrices, Si = e 2lSl , in these two channels. 

Note that a bound state has the most profound influence on phase shifts in the limit 
Eb | 0. In this limit 

+ + (12) 
and the phase-shift flip occurs (cf. |18|). On contrary, in the limit | — oo, 

Si^^(\l\-\l + a\). (13) 

3 The Krein-Friedel formula and the DOS 

The DOS in the presence of the AB potential is defined to be 

p a {E) = ImTr G Q (x, x, E + ie), (14) 

7T 

where G a {x,y, E + ie) the resolvent (the Green function) of H. The integrated density 
of states N a (E) is then as usual given by 

N a (E) = I* p a (E')dE'. (15) 

J — oo 

The DOS in two dimensions when all interactions are switched-off is Po(E) = (m/2nh 2 )V, 
with V = J d 2 r being the (infinite) volume. To calculate the change of the integrated 
density of states (IDOS) in the whole space we shall make use of the Krein-Friedel formula 
Q that gives the change AN a (E) of the IDOS induced by the the presence of a scatterer 
of a finite range directly by summing over phase shifts, 

AN a (E) = N a (E) — N (E) = -Y) = (2m)- 1 In det S, (16) 



with S the total on-shell S-matrix. The fact that phase shifts can be rather easily calcu- 
lated without any care of the proper normalization of wave functions greatly facilitates 
the calculation. Moreover, by means of the Krein-Friedel formula it is rather easy to 
calculate the change of the IDOS for all possible self-adjoint extensions. In the case of the 
long-ranged AB potential we have found that the Krein-Friedel formula when combined 
with the ^-function regularization can still be used despite the fact that phase shifts (|J) 
are in general singular [EL U\. In the absence of bound states, 



OO 



lndetS= 2i$i = ™ E ~ \ l + a \) 



Z7T 



DO 



2E^-E(^r-E('-^ 

. 1=1 1=0 1=1 



s=-l 



= iir [2( R {s) - Ch(s, rj) - ( H (s, 1 - r ? )]| s= _ 1 = -iwq[\ - rj), (17) 

where (r and (h are the Riemann and the Hurwitz ^-function. Thus, using (|16|), the 
change of the DOS is 

Ap a (E) = Pa (E) - p (E) = - 1 - 77(1 - r?) 5(E), (18) 

and Ap a (E) is only the function of a distance from the nearest integer. 

In the presence of bound states, the contribution of scattering states to AN a for E > 

is 

a , T / 1 / i 1 / sin(?77r) \ 

AN a (E) = - —77(1 - r]) + - arctan ' 1 



'/ 



2' V " IT V COS (^) - (\ E -n\/E) 

1 / sin(?77r) \ 

- ^^Uw+dB-.-iw-'J' (19) 

where E_ n and E_ n _i are the binding energies in I = —n and I = —n — 1 channels. By 
repeating the same calculation for a < one finds that AN a (E) is a symmetric function 
of a §, 

AN. Q (E) = AN a (E). (20) 



3.1 The resonance 

Note that for < r\ < 1/2 the resonance appears at 



The phase shift 5- n (E) ([0]) changes by ir in the direction of increasing energy and the 
integrated density of states ([H]) has a sharp increase by one. The profile of the resonance 
[the argument of arctan in (|TT])] is given by 



E r > - E% s E^ - El 



v ' 



22 



where T = i±^ es tan 7/7r is the width of the resonance. Note that its profile (|22|) is not of 



the Breit-Wigner form (see Ref. |19 |, & 145), 

r 



(23) 



E — E 

For 1/2 < 7] < 1 the resonance is shifted to the I = —n — 1 channel. 77 = 1/2 is a special 
point since resonances occur in both channels at infinity. Therefore the contribution of 
the arctan terms in ([19]) does not vanish as E — > 00, but instead gives —1. 

4 Regular izat ion, R — ► limit, and the interpretation 
of self-adjoint extensions 

Different self-adjoint extensions correspond to different physics inside the flux tube (see 
an example in Ref. P0|, p. 144). To identify the physics which underlines them we have 



considered the situation when the AB potential is regularized by a uniform magnetic field 
B within the radius R and satisfies the constraint 



B(t) d 2 r = $ = const. (24) 

n 



One finds that in the absence of the magnetic moment (g m = 0) or any other attractive 
interaction with the interior of the flux tube the matching equation for the exterior and 
interior solutions in the l-th channel is (see for example), 



_ , ,,, . l*l + j + 1 + (* 2 A*) iFi (ziz + ( x2 / 4q )' l z ' + 2 ' a ) 

X K M (x) a + ^ + a + 1 1 F 1 (^±i + (xV4a),|/| + !,«)' { ] 



Here i-F\(a, 6, c) is the Kummer hypergeometric function |22|, and x\ = K[R ^ 0. However, 
since the l.h.s. decreases from — \l + a| to — oo as x — > oo and the r.h.s. is always greater 
than — I + \a\ one finds that Eq. ( f25|) does not have a solution unless it is an attractive 
potential V(r) inside the flux tube, 



where V(r) = otherwise. Here, c(R) = 2(1 + e(R)), s(R) > 0, and e(R) -> as R -> 
2~T, p3j . This amounts to changing x 2 /2a to x 2 /2a — c/2 on the r.h.s. of (pop. Note that 



in the limit -R — > 

V(r)\ r<R ^---S(r). (27) 
m r 

The attractive potential can be either put in by hand or, if the Pauli Hamiltonian is 
considered, as arising from the magnetic moment coupling of the electrons with spin 
opposite to the direction of the magnetic field B. In the latter case the critical potential 
corresponds to the case when gyromagnetic ratio g m = 2 [e(R) = (g m — 2)/2 = 0] (cf. 
Eq. (H). Then the matching equation (|25f) has a solution in x = for / = —n. Since the 
magnetic field is not singular any more the Aharonov-Casher theorem || applies. It is 
known that there are ]a[— 1 zero modes in this case, }a[ the nearest integer larger than 
or equal to a. If a is an integer then one has exactly n — 1 zero modes, if not, their 
number is n = [a] [Q. The result only depends on the total flux a and not on a particular 
distribution of a magnetic field B. 

Whenever g m > 2 (and hence e > 0) or g m = 2 with an attractive potential V(r) = 
—e/R 2 , e > arbitrary small, the bound states may occur in the spectrum in the channels 
I < 0. They correspond to solutions xi > of In other words the coupling with the 
interior of the flux tube becomes sufficiently strong for the particle to be confined on the 
cyclotron orbit inside it. Note that the wave function (0) of bound state decays expo- 
nentially outside the flux tube. In contrast to the zero modes their number does depend 
on a particular distribution of the magnetic field B. Using three different regularizations, 
uniform, regular, and cylindrical one finds that the number of bound states is less than 
or equals to g 0, |4j 

# 6 = 1 + n + [a(g m - 2)/4] + [a(g m + 2)/4 - n] , (28) 



with [.] as above. Note that the number of bound states is generally higher than the num- 
ber of zero modes 0. The bound is saturated || 0] if the cylindrical shell regularization 
frSfl of the AB potential is used. The physical origin of this difference can be understood 
in a simple way. In the latter case the energy Eb of magnetic field is infinite for any 
R and in this sense the magnetic field inside the flux tube is much stronger than, for 
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example, in the homogeneous field regularization when Eb, 

E b = nB 2 R 2 /2 = $ 2 /(27ri? 2 ), (29) 
stays finite for any nonzero R. 

5 The Hall effect 

The S matrix, s a (ip), in the AB potential was calculated according to 

1 CO 

S M = z~ E e 2l5l+Uv - (so) 

l=— oo 

We have found that either in the absence or presence of bound states, 

8-\ a \(<p) = s\ a \(-<p), (31) 

under the transformation a — *> —a 0. 

In the presence of bound states, the differential scattering cross section for <p ^ was 
found to be @ 

4 sin(7rcM 

+ t ■ ( /9 x [sin A_„cos (A_ n - 7ra + (p/2) + sin A_„_i cos (A_ n _i + na - (p/2)] , (32) 

K Sllll I Z I 

where 

( dcr°\ 1 sin 2 (vra) 

( fc ^) = o i • 27 /o\ ( 33 ) 



dcp J ' 2irk sin ((p/2) 



is the differential scattering cross section in the absence of bound states |L6| . The peri- 
odicity of the differential cross section with respect to the substitution a — > a ± 1 then 
follows from Eq. fl3"2|). Note that in the presence of bound states, the differential cross 
section becomes asymmetric with regard to (p —> —tp (what is equivalent, with regard to 
a — > —a). The origin of the asymmetry is easy to understand since bound states for 
a > are only formed in channels with I < 0. 

Now, if one considers a random distribution of flux tubes with the density n v , then 
the Hall effect is induced H]. The Hall resistivity p xy can be calculated in the dilute 
vortex limit, i. e., when the multiple-scattering effects are ignored. The quantity that 
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measures the fraction of the electrons moving in a transverse direction is sin(yj) da(kp, ip). 
Therefore, if the density of vortices is n v , the Hall current in the dilute vortex limit is 
proportional to 

d(p sin (p — (kF, <p)- (34) 
-7T dip 

By inverting the conductivity tensor one finds that the Hall resistivity, p xy , is 

k F r dip , da 



Pxy 



o k F r dip . da 

Ph— 7r smi ?l~( kF >W> ( 35 ) 
a J-ttZti dip 



which was also obtained by Nielsen and Hedegaard |25[. Substituting ( p^ ) to (^) then 
gives 

An he 

Pxy = — — — n sin(7ra) [sin A_ n cos(A_ n — ira) + sin A_„_! cos(A_„_! + 7ra)l . (36) 
n e e 1 

The result shows that the Hall resistivity is proportional to the density of vortices and 
depends on their vorticity via trigonometrical functions. In particular, as a self-consistency 
check, the Hall resistivity (|3~6"D vanishes for a an integer. The Hall resistivity also vanishes 
whenever A_ n = — A_ n _i modulo it. 

6 Discussion of the results 

The single-particle density of states p a induced by the AB potential was calculated. It was 
shown that p a is a symmetric and periodic function of the flux. Existence of the resonance 
in the AB potential was proven and the phase-shift flip was discussed. In addition to the 
flux, the number of bound states for a nonsingular flux tube was also shown to depend 
on the energy of the magnetic field. The Hall resistivity in the dilute vortex limit was 
calculated. 



Our results are not only of academic but also of practical interest |13| thanks to the 



recent developments in the fabrication of microstructures and in mesoscopic physics (see 



[ 26 1 for a recent review). In particular, challenging is observation of the resonance and the 
Hall resistivity. They occur only in the case if a bound state is present or, in the latter case, 
if the phase-shift flip occurs. It will be interesting to consider an application of our results 
in the set-up (see Fig. 0) proposed by Rammer and Shelankov JT2|, especially because 
recent measurements on yttrium-barium-copper oxide (YBCO) delta rings with three 
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grain-boundary Josephson junction | 2"?fl , reported the observation of vortices that curry a 



flux a = 1/4 which is smaller than the standard flux quantum hc/2e (corresponding to 
a = 1/2) in the superconductor. Therefore, when the high-T c YBCO film is used as a 
gate on top of the heterostructure containing 2DEG, the resonance is at some finite energy 
and could in principle be observed. The same experimental set-up is also promising for 
detecting the Hall effect for the random distribution of vortices. 

Our results have been presented in the mathematical language of self-adjoint exten- 
sions. A self-adjoint extension is actually the R — > limit, where R is the radius of a flux 
tube. Experimentally, infinitely thin means nothing but that the radius of the flux tube is 
negligibly small when compared to any other length, such as a wavelength of particles, in 
the system. Therefore, this is the regime in which our results can be applied. Parameters 
A_ n _ n _i of self-adjoint extensions are then determined by bound state energies in the 
/ = — n and I = —n — 1 channels. 

I should like to thank A. Comtet, Y. Georgelin, S. Ouvry, and J. Stern for many useful 
and stimulating discussions. 
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